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SUBSTRUCTURE COUPLING FOR DYNAMIC 
ANALYSIS AND TESTING 

Roy R. Craig, Jr.* and Ching-Jone Changt 
The University of Texas at Austin 


SUMMARY 


Substructure coupling methods may be classified as fixed-interface 
methods, free-interface methods, or hybrid methods, depending on whether the 
modes used are determined with substructure interfaces fully restrained, 
free of restraint, or partially restrained. 

In the present report three methods are presented for reducing the 
number of interface coordinates required by fixed interface methods. These 
are: Ritz reduction, Guyan reduction, and modal reduction. Example solu- 

tions are given to illustrate these three methods, and procedures for estima- 
ting convergence properties of the methods are described. 

Also presented in the report is a study of a class of free-interface 
methods. The methods described appear to be significantly more accurate than 
the free-interface methods of Hou and Goldman and to be particularly useful 
when experimental verification of the analytical model is required. The sub- 
structure equations and coupling equations are developed, and example prob- 
lems are presented which compare the accuracy of the present methods with 
other free-interface methods and with fixed-interface methods. 

♦Associate Professor, Aerospace Engineering and Engineering Mechanics 
tGraduate Student, Engineering Mechanics 



INTRODUCTION 


Substructure coupling methods may be employed when (1) the number of 
coordinates in a dynamic analysis needs to be reduced, or (2) when analysis 
or testing of a structure needs to be done separately on different portions 
of the structure. Substructure coupling methods generally employ coordinate 
transformations where the generalized coordinates of substructures are defined 
by Ritz-type modes. In most cases some, or all, of the generalized coordi- 
nates are defined by a set of substructure free vibration modes. The sub- 
structure coupling methods may be classified as fixed-interface methods, free- 
interface methods, or hybrid methods depending on whether the modes are deter- 
mined with the substructure interfaces fully restrained, free of restraint, 
or partially restrained. 

The primary focus of the research described in the present report is 
on methods for reducing the number of interface coordinates required by 
fixed-interface substructure coupling methods. The fixed-interface method 
developed by Hurty (ref. 1) and refined by Craig and Bampton (ref. 2) requires 
that all of the original interface displacement coordinates be retained in 
the final coupled structure. Benfield, Bodley, and Morosow (ref. 3) conducted 
a study comparing the accuracy of several substructure coupling procedures. 
Figure 1 shows some of the results obtained in reference 3. Although the 
various methods were only applied to one structure, the results of this 
limited study indicated that the fixed interface method of Hurty and the 
equivalent Craig-Bampton method generally produced the most accurate solutions 
for a given number of degrees of freedom. Slightly better accuracy was 
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achieved by methods introduced by Benfield and Hruda (ref. 4) but these 
latter methods suffer the disadvantage that the substructure modes of one 
substructure are not independent of the other substructures. Benfield, et. al. 
(ref. 3) suggest that a drawback to the Hurty method (and Craig-Bampton 
method) is the requirement that all interface displacement coordinates, i.e., 
constraint modes, must be retained iri the final coupled structure. It is 
suggested that this might restrict the number of component, or substructure, 
normal modes retained in the coupled structure. It is the purpose of this 
report, therefore, to describe three methods for reducing the number of inter- 
face coordinates in the final system model. Reference 5 presents additional 
information concerning these three methods. 

A secondary purpose of this report is concerned with improving the 
accuracy obtainable from free-interface substructure coupling methods. Figure 
1, taken from reference 3, indicates that the two free-interface methods, 

Hou's method and Goldman's method, produce frequencies much less accurate 
than those produced by fixed-interface methods for a comparable number of 
system degrees of freedom. 

Subsequent to the appearance of reference 3, Rubin (ref. 6) presented a 
paper describing a class of free-interface methods similar to ones described 
previously by MacNeal (ref. 7) and Klosterman (ref. 8). Although Rubin went 
into considerable detail to describe how to represent individual unconstrained 
(i.e., possessing rigid-body freedom) substructures, he did not fully develop 
the substructure coupling procedures or present any examples of multi -sub- 
structure systems. 

Rubin (ref. 6) stressed the fact that, when experimental verification 
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of analytical models is required, the class of free-interface substructure 
methods which he presented are much more attractive than are fixed-interface 
substructure methods. 

Since the methods presented by Rubin appeared to have the potential for 
producing substructure coupling results of acceptable accuracy and appeared 
to be compatible with state-of-the-art dynamic testing procedures, it was 
decided that the initial effort along the lines of improving free-interface 
methods should be directed toward further development of the methods suggested 
by Rubin. Rubin's own work (ref. 6) treats only unconstrained substructures. 
The only example presented is for a single substructure represented by a con- 
tinuous (i.e., partial differential equation) model. In the present report 
the theory is presented first for constrained substructures, in order that 
the concepts may be presented first in the simplest form possible. Also, in 
the present report only discrete models (i.e., matrix ordinary differential 
equations) are employed, since these are the ones employed to model practical 
structures. Finally, the present report presents detailed equations for 
coupling substructures and presents example problems in which results obtained 
by Rubin-type methods are compared with results obtained by other methods 
(Hou, Craig-Bampton) . 
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FIXED INTERFACE METHODS WITH REDUCTION 
OF INTERFACE COORDINATES 

In this part of the report three methods are derived for reducing the 
number of interface coordinates required for fixed-interface substructure 
methods. Examples are also presented to assess the accuracy of the derived 
methods . 


Substructure Equations 

For a typical substructure, the equations of motion, assuming linear 
viscous damping, may be written 

mx + cx + kx = f (1 ) 

For undamped, free vibration, the equations reduce to 

mx + kx = f (2) 

where f contains forces at junction (interface) coordinates in the coupled 
structure. 

In most substructure coupling methods the original substructure physi- 
cal coordinates are expressed, after the manner of Ritz, in terms of a set 
of substructure generalized coordinates. In the method of Craig and Bampton 
(ref. 2) these substructure generalized coordinates consist of a set of junc- 
tion, or interface, coordinates and a set of normal mode coordinates. These 
normal modes are obtained with all junction coordinates fully restrained. 
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Equation (2) may be partitioned to give 
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where x i contains the "interior", or non-interface, physical coordinates and 

x- contains the "junction", or interface coordinates. Figure 2 shows interior 
J 

nodes and junction nodes of two adjacent substructures. 

Since the normal modes are to be obtained with x^ = 0, Eq. (3) reduces 
to 


The mode shapes, 4> n , can be assumed to be normalized so that 


<fc m d> = 1 

v n ^n 


*n M n = X n 


(5) 

(6) 


and the mode shapes collected to form the respective columns of a modal 
matrix The substructure "frequency" has been designated x to distin- 

guish it from the system frequency w. Later, the modal matrix will be parti- 
tioned so that some of the modes are kept in the final analysis while other 
modes are reduced out. 

If there were no loads (including inertia loads) at the interior 
coordinates of the substructure, the static displacement would be given by 
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( 7 ) 
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Displacements x- of the junction coordinates would produce displacements of 

J 

the interior coordinates as given by the upper partition of Eq. (7), namely 
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The components of the transformation matrix can be obtained by solving 

’ J 

the equation 




( 10 ) 


A substructure coordinate transformation may be defined now by 


f 


— 

— 



X . 


. 

<&. • 


p. 

1 


11 

1J 


K 1 

< 

► — 



J 

> 

X . 


0 

I 


p . 

J 






k J 


— 

_ 


W J 


or 


x = $ p 
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( 12 ) 


The equations of motion may be transformed to the new substructure generalized 
coordinate system as follows: 


mp + cp + kp = g 


(13) 
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where 


m = $> m $ 

c = $ T c $ 

k = $ T k $ 

T , 

g = <j> f 


( 14 ) 


The transformed matrices may be written in more explicit form through 
the use of Eqs. (11) and (14). Thus, 





The null partitions of R result from the transformation imposed by Eq. (10). 
The remaining submatrices are given by 
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Coupled System Equations 


In discussing the coupled system, we will consider two substructures, 
"A" and "B". Extension to more substructures is straightforward. The total 
set of system generalized coordinates is 


q 



( 17 ) 


where is a collection of all substructure (interior) normal mode coordi- 
nates and q. is a set of all junction coordinates. Although it is common 
practice at this point to reduce the number of coordinates by including only 
selected modes from each substructure, we will assume that all modes are 
retained for the time being and will treat the reduction of coordinates 
question in a later section. 

The substructure coordinates are related to the system coordinates by 
"label" submatrices, which simply identify substructure coordinates with the 
appropriate system coordinates. This is a formal expression of the "direct 
stiffness" procedure for assembling system matrices from element matrices. 
Thus, 
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or 

P = L q . (19) 

The system matrices are obtained by direct stiffness assembly procedures 
using the labelling information of Eq. (18). The resulting system equations 
may be written 

Mq + Cq + Kq = F . (20) 

In partitioned form the equations of motion for free vibration of the 
undamped coupled system are 



For the undamped free vibration problem we can take 

q(t) = Q cos (cot) (22) 
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where Q is a vector of amplitudes only. The resulting eigenvalue problem 
has the form 



where the substructure normal mode coordinates have been combined into a 
single vector cos (wt). 

In the following sections we will discuss the effect of reducing out 
some of the substructure normal mode coordinates and/or some of the junction 
coordinates. 




Reduction of Coordinates 


One of the primary reasons for adopting a substructure method for solv- 
ing structural dynamics problems is to reduce the number of degrees of freedom 
which must be employed in the dynamic analysis. The truncation of substruc- 
ture normal modes in a fixed-interface substructure coupling analysis has been 
studied by Hurty (ref. 9), Benfield, Bodley, and Morosow (ref. 3), and others. 
For completeness, this procedure is described below. The reduction of junc- 
tion (interface) coordinates has apparently not been considered previously. 
Three possible methods for reducing junction coordinates are introduced in 
succeeding sections. These are referred to as modal reduction, Guyan reduc- 
tion, and Ritz reduction. 

Reduction of substructure normal mode coordinates . - For considering 
reduction of normal mode coordinates we can partition Q i in Eq. (23), giving 
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where coordinates Q g will be retained and coordinates Q b will be reduced out. 
The second row-partition of Eq. (24) may be written 

(K bb - “ 2 M bb> «b = “ 2 M bj (25) 

or 
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(26) 


“ 2 <’ - “ 2 K bb H bb>~ ’ K bJ % Oj 


Consider the eigenvalue problem 


(K bb ■ x V Q b = 0 < 27 > 

i.e., the eigenvalue problem associated with the substructure modes to be 

reduced out. The frequencies of these substructure modes are called A to 

n 

distinguish them from frequencies o> of the coupled system. Recall that the 

modes are normalized so that M bb = 1 and is a diagonal matrix containing 
2 

A n 's. Thus, Eq. (26) takes the form 


(w/A ) 

Qb 1 - (o)/A n ) 2 Mbj Qj 


Let A ] be the lowest frequency among the eigenvalues of Eq. (27), i.e., 
the lowest frequency of a mode which is to be reduced out. If the range of 
system frequencies, co, of interest is such that 


w/A-j << 1 


the zeroth-order approximation to Q. is 


Q b - 0 


i.e., the modes associated with Q b can just be eliminated. The reduced 
eigenvalue problem is simply 
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(31) 
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An iterative procedure for improving frequencies and mode shapes of the 
system could be based on using the zero-order approximate values of u in 
Eq. (28). Kuhar and Stahle (ref. 10) proposed a "dynamic transformation" 
procedure which essentially consists of defining in terms of Qj through 
Eq. (28) using a pre-selected transformation frequency, p, instead of the 
(unknown) system frequency o>. However, involves mode shapes of the 
deleted modes, and hence this procedure requires that such modes be obtained 
in the original substructure eigensolutions. 

Examples will be presented in a later section to compare the effective- 
ness of reducing out substructure normal modes with the effectiveness of 
reducing out junction coordinates. 

Modal reduction of junction coordinates . - This section will introduce 
the technique of modal reduction of junction coordinates. Consider the 
eigenvalue problem that would be obtained from Eq. (23) by constraining 
the substructure normal mode coordinates and allowing motion of junction 
coordinates, i.e., 

< K jj - - 0 < 32 > 


It may be noted that this corresponds to using Guyan reduction to reduce out 
all interior coordinates. The frequencies of these "junction modes" will be 
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called <5 to distinguish them from system frequencies, w, and substructure 


frequencies, X. Let the junction modes form a modal matrix 4'. , which has 
the form 
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where the modes 4^ c will be kept and modes 4^ 
that the modes are normalized with respect to 
Now let 
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M.. in the usual fashion. 
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or 

Q = 4- Q 


(36) 


Transformations of the form 

K = y T K v 


(37) 


transform Eq. (23) to the form 
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(38) 


The zeros in the c-d coupling terms are due to orthogonality of the junction 
modes. The submatrices in Eg. (38) are obtained from the equations 
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Let Eq. (38) be written as two equations, namely 
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Since, according to Eqs. (39), and are diagonal, Eq. (41) may be 
written 
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where 


“dd = P w nJ 

is a diagonal matrix whose elements are given by 
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When Eq. (42) is substituted into Eq. (40) the resulting equation may be 
written 
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where 


AM.. = M. , W,, M , . = Jw M. . ip ip M-. 
n id dd di L n ij v n ji 


(46) 


in which i|> n is the junction mode corresponding to 6^ . 

The derivation of "modal participation factors," which provide an 
estimate of the contribution of each junction mode to the solution of Eq. (45) 
proceeds along lines similar to the procedure employed by Hurty in reference 
9 for discussing truncation of substructure normal modes. 

First a "base problem" is defined by setting aM^. = 0 in Eq. (45) or, 
equivalently, setting = 0 in Eq. (40), giving 
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The exact solution of Eq. (38) and the solution of Eq. (47) are related by 
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When Eqs. (47), (48), and (49) are substituted into Eq. (45) and "higher 
order terms" in Au are neglected, there results 
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If Eq. (51) is premultiplied by [QT q qJ q ], aQ is expanded in terms of the 
eigenvectors, Q q , and proper account is taken of Eq. (47) and orthogonality, 
there results 
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If the solution vectors of the base problem are normalized relative to the 
mass matrix of the base problem the denominator of Eq. (52) becomes unity, 
so the eigenvalue error is given by 




AMi i ( 
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But, since AM.... can be written as a sum, as shown in Eq. (46), Eq. (53) 
may be written 
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From Eq. (54) it is clear that modal participation factors can be defined by 
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(55) 


Two factors are of interest in determining the value of the modal 
participation factor: (a) the frequency ratio, , and (b) the mass 

coupling. If u Q /6 n « 1 the nth junction mode should have little contribu- 
tion to the frequency of a given system mode. On the other hand, Eq. (55) 
shows that the coupling of the interior portion of the base problem mode, 
i.e., Q i0 , with the nth junction mode through the original coupling mass 
matrix M^. is also important. Hence, the modal participation factor given by 
Eq. (55) should, where possible, be used as a basis for selecting junction 
modes to be retained. 

It is clear from Eq. (55) that the concept of modal participation factor 
has meaning only in reference to choosing modes which would lead to a maximum 
improvement upon a given base solution. The question arises, then, as to 
what should constitute the base problem, i.e., what junction modes belong to 
^jc 1n (34). It would seem that the only rational approach to this prob- 
lem would be to select a specified number of junction modes using junction 
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frequency alone as the criterion, i.e., order the junction modes according to 
frequency and use the modes from the lower end of the frequency spectrum. 

If there are any rigid-body junction modes, these must be included in the base 

problem. 

In a later section, an example will be presented to illustrate the 
calculation of modal participation factors. 

Guyan reduction of junction coordinates . - In reference 11, Anderson, 
et. al . , applied the method of Guyan reduction to a plate vibration problem 
and showed that acceptable accuracy of frequencies was maintained even when 
a large percentage of the coordinates were reduced out. In this section, a 
similar Guyan reduction will be applied. However, in this case only selected 
junction coordinates will be reduced out. 

Consider again Eq. (23) and partition Qj into 



where the coordinates are to be reduced out by a Guyan-type reduction. 
Thus, Eq. (23) may be written 
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Guyan reduction is based on neglecting inertia terms in the third row-parti- 
tion of Eq. (57) to obtain 


where 
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Although there is no established procedure for selecting coordinates to 

be retained and coordinates to be deleted, i.e., for ordering the Q. coordi- 
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nates in Eq. (56), a possible procedure for selecting the coordinates to be 

retained might be based on magnitudes of the quotient of terms on the diagonal 

of K.. divided by corresponding terms on the diagonal of M--, with the coordi - 
JJ J J 

nates having the smaller k/m values being retained. 

The junction coordinate reduction equation, Eq. (59), has the same form 
as a substructure equation, Eq. (8). The system displacement transformation, 
based on Eq. (58), is 
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Note that the coordinates Q f are reduced out, i.e., approximated by Eq. (58). 
They are not simply eliminated as were Q b , Eq. (30), and Q d , Eq. (47). 
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Transformation of Eq. (57) to the reduced set of coordinates Q 


defined by Eq. (60) gives 
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Note the similarity of the above expressions to those in Eq. (16). 

An example of the use of Guyan reduction will be presented in the 
section on example problems. 


Ritz reduction of junction coordinates . - The modal reduction and Guyan 
reduction of junction coordinates which were presented in the preceeding 
sections are both special cases of Ritz representation of junction coordinates. 
There are cases where the nature of the problem might suggest a direct 
approach of Ritz at either the substructure level or the coupled system level. 
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We will employ the latter approach, beginning with Eq. (23), repeated here 
for convenience. 
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Let the junction coordinates be expanded in terms of a set of Ritz 
vectors according to the equation 
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where notation similar to that of the preceding sections is employed. The 

matrix 4^ is rectangular, (nj > n^), since the number of coordinates is to 

be reduced by this Ritz expansion. The columns of 4'. represent the displace- 
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As before, this can be represented by 

Q = ^ Q 

and the stiffness matrix becomes 
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Then the system equation, Eq. (23), is transformed into 
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where K... and are diagonal. 

The use of Eq. (67) for reduction of junction coordinates will be 
illustrated in the section on example problems. 

It is possible to write expressions for error analysis, or for modal 
participation factor analysis, for the general case of Ritz representation 
of junction coordinates. First, Eq. (63) will be expanded by assuming that 
one or more Ritz vectors (or modes) are to be added in Eq. (63), with the 
set Yjg forming the "base problem" and being the added vectors. Thus, 
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This gives a base problem similar to Eq. (67), namely 
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and an "augmented problem" of the form 
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where, in addition to terms defined by Eq. (68), 
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Equation (71) can be written 
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There might be instances where the matrix to be inverted in Eq. (74) is 
singular, but this would simply mean that the Ritz vector causing the singu- 
larity is very important at the given frequency w, since this would imply a 
very large Q h - The matrix inversion in Eq. (74) leads to the definition of 
w n in Eq. (44) when junction modes are used as the Ritz vectors. 

Equations (73) and (74) may be combined to give 
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To simplify notation, let Eq. (75) be written 
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where 


where 
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Let the solution of the "augmented problem," Eq. (71), and the solution 
of the "base problem," Eq. (70), be related by 
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Q 0 + Q 


(79) 


and assume that the eigenvectors of the base problem are normalized on the 
mass matrix of that problem. 

For the case of modal junction coordinates it is easy to simplify the 

expressions given in Eq. (77) above. Certain assumptions on the expansion of 

2 

Eq. (77) in powers of go allow the same to be done here. For purposes of 
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error analysis, i.e., determining modal participation factors, let it be 
assumed that the S matrices in Eq. (76) can be approximated by S(co Q ) giving 
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Equation (80) is analogous to Eq. (51). Equations (70) and (79) were employed 
in converting Eq. (76) to Eq. (80). If Eq. (80) is premultiplied by qJ , aQ 
is expanded in terms of the eigenvectors of the base problem, and proper 
account is taken of Eq. (70) and orthonormality of the eigenvectors of the 
base problem, then Eq. (80) reduces to 
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Equation (81) reduces to Eq. (53) for the case of modal junction coordinates. 
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If a single Ritz vector is used in ¥- h in Eqs. (72), then Eq. (81) becomes a 
modal participation factor" determining the approximate error produced by 
omission of that vector from the base problem, or, conversely, the approximate 
improvement of the base problem solution which could be achieved by augmenting 
the solution by that Ritz vector. 
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Examples 


A limited number of example problems have been solved to illustrate 
methods for reducing junction coordinates. Three structures were employed in 
the studies: the two-substructure truss of reference 3, a nine-substructure 

ring, and a two-substructure grid. The first two structures possess rigid- 
body freedom; the third is a cantilever structure. 

Dimensions have not been given in the example problems of this report 
since the purpose of the report is to compare relative accuracies of various 
methods of solving given problems. 

Truss examples - modal reduction and Guyan reduction of junction 
coordinates . - The truss shown in Figure 3 was used by Benfield, et. al. 

(ref. 3) in their comparison study of substructure coupling methods. Results 
will be presented for both modal reduction of junction coordinates and Guyan 
reduction of junction coordinates of this truss. 

Table 1 gives results for frequencies of the above truss obtained 
by the use of modal reduction of junction coordinates. The column headed 
"Exact" gives results for the coupled structure with all of its original 
degrees of freedom retained. The remaining results were obtained using 
substructure coupling with 5 normal modes retained for each substructure, 
i.e., = Ng = 5. The columns headed with Nj = 6 give results based on the 

original Craig-Bampton method, since N- = 6 is the total number of junction 

J 

coordinates. Since the truss as a whole is free to execute 3 degrees of 
rigid-body motion, the columns headed N- = 3 give results based on assuming 

J 

the junction to be rigid, i.e., the 6 junction coordinates must move as a 
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rigid body. Nj = 3 is the smallest number of junction coordinates that may 
be retained, since it is essential to allow at least rigid-body motion of the 
junction. The columns headed N • = 5 and N • = 4 give results based on the 

J J 

modal reduction procedure. Note, for example, the frequencies for elastic 
mode 3. The accuracy of the frequency is not seriously degraded by reduction 
of junction coordinates until the last stage (N- = 3), but mode 3 does require 

J 

more than rigid-body behavior from the junction coordinates. The choice of 
junction modes to be retained was based on junction mode frequency alone. 

In Table 1 it may also be noted that the errors in some modes, e.g., 
elastic modes 7 and 8, are due to truncation of substructure normal modes. 

The large errors persist even when all six junction coordinates are retained. 

Table 2 presents data on substructure normal mode frequencies and 
junction frequencies. It may be noted that the junction frequencies are the 
same order of magnitude as the higher-frequency substructure normal modes 
(e.g., 5th through 10th modes). 

It was previously noted that junction mode 4 (S 2 = 0.0182248) makes a 
significant contribution to system elastic mode 3 (w 2 = 0.0030508). Table 3 
presents modal participation factor data for this mode so that an assessment 
may be made of the relative importance of frequency ratio, a) 0 /6 n , and the mass 
coupling term as they appear in Eq. (55). 

The results presented in Table 3 illustrate the fact that both the 
junction frequency ratio and the mass coupling factor are important in deter- 
mining the modal participation factor. In the present example it is the low 
value of the mass coupling factor for junction mode 5 that renders it unimpor- 
tant in the solution for this particular system mode. The estimated error 
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Table 1. Frequencies of 2-Component Benfield Truss with Modal Interface Reductions 
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Junction Modes in MPF 


n = 4 

n = 5 

n = 6 

<»6>o 

r2 / \ 2 

6 “ [Ur) 

n v 6 q 

0.2268 

0.1098 

0.0856 

*] 0 M ij *n 

0.7121 

5.8 x 10' 14 

0.1450 

p - A n ( “6> 2 

" " K>o 

I 

-0.1150 

o 7 to-28 

-3.7 x 10 

-1.8 x 10' 3 


Table 3. Modal Participation Factor Data for Mode 6 
(Elastic Mode 3) of Benfield Truss - 
Base Solution: Nj = 3, = Ng = 5; 

UJ 3 = 0.0033692 ; (u>,) 2 ^ = 0.0030508 

6o v 6 'exact 
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due to omission of junction mode 4 from the base solution is 11.50%; the 

actual error was 10.44% (see Table 1). 

The Benfield truss was also used to illustrate the Guyan reduction 

method. The numbering of the six junction coordinates is shown on Fig. 3. 

The solution obtained by retaining coordinates 1, 2, and 3 is the same as the 

N. = 3 solution by the modal reduction method, since this case consists of 
J 

modelling the interface as a rigid-body. Guyan reduction was employed by 

first reducing out coordinate 6 (leaving Nj = 5), and then reducing out both 

5 and 6 (N. = 4). Had other choices been made of coordinates to reduce out, 

J 

the results would probably have been different from those shown in Table 4. 

The results of Tables 1 and 4 are summarized in Fig. 4, which shows a 
comparison of the number of modes having an accuracy of 0.5% or better. 
Although modal reduction of interface coordinates yielded slightly better 
results than Guyan reduction, the present example is too small for any general 
conclusions to be drawn from it. 

Grid examples - modal reduction of substructure normal mode coordinates 
and junction coordinates . - Figure 5 shows a cantilever grid structure used 
to illustrate further the modal reduction of junction coordinates. A total 
of 12 junction coordinates couple the two substructures together. 

This grid structure was used to provide an example having more junction 
coordinates than the truss previously discussed and having no rigid-body 
freedoms. Studies were made to compare the effect of reducing junction mode 
coordinates with the effect of reducing substructure normal mode coordinates. 
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Table 5 gives the frequencies obtained using various combinations of 
substructure normal modes and junction modes. In all cases, the modes to be 
retained were selected on the basis of frequency alone. Figure 6 shows a 
comparison of the number of modes having an accuracy of 0.5% or better. 

Figure 6 seems to indicate that modal reduction of junction coordinates and 
reduction of substructure normal mode coordinates produce similar decreases 
in the number of modes accurate to within 0.5%. Again, the example problem 
is too small to permit generalizations concerning the effect of modal reduc- 
tion of junction coordinates, but the results shown in Fig. 6 suggest that, 
within limits, modal reduction of junction coordinates is feasible. 

Since there are only 4 translational coordinates along the junction, it 
might appear that 4 junction mode coordinates would be adequate to represent 
junction behavior in a good many, if not all, modes. Case 3 (N A = 12, 

N B = 8 ’ N j = 4 ) indicates tha t this is not so. Note that Mode 9 has a large 
frequency error (12.03%) in Case 3 but only 0.3 % error in Case 2. It was 
discovered that the 5th junction mode, i.e., the last one eliminated in going 
from Case 2 (Nj = 8) to Case 3 (Nj = 4) is predominantly torsion of the 
junction line, and that this mode contributes significantly to both 9th and 
10th system elastic modes. 

The above discussion serves to caution against reducing out junction 
modes using a frequency criterion alone. Modal participation factors, as 
previously described, should, if possible, be used in choosing which modes to 
retain. Also, whereas previous studies of Guyan reduction of system coordi- 
nates have indicated that most, if not all, rotation coordinates could be 
reduced out, it is clear from the present example that junction modes in 
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Table 5. Frequencies of Cantilever Grid Structure 
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Table 5 (cont.). Frequencies of Cantilever Grid Structure 




which "twist" of the junction line is significant may be important to the 
accuracy of system frequencies obtained through the use of modal reduction of 
junction coordinates. 

Ring examples - Ritz reduction of junction coordinates . - Figure 7 shows 
a ring structure which was used to provide examples of Ritz reduction of junc- 
tion coordinates. The ring was divided into nine substructures, with each 
substructure having eight elements. Since each element subtends an arc of only 
five degrees, a straight beam element was used, with the translational coordi- 
nates transformed to the radial and circumferential orientation at the nodes. 
Although they participated very little in the solution for the first ten 
system modes, two substructure normal modes were included for each substruc- 
ture in generating the substructure coordinates in Eq. (17). 

Two cases are presented: (1) a uniform ring, which possesses symmetry 

properties and for which an analytical solution is available, and (2) a non- 
axisymmetric stiffened ring. In each case solutions employing a Ritz approxi- 
mation of junction coordinates based on Eqs. (63) and (67) are compared with 
solutions employing all 27 junction coordinates plus 18 substructure normal 
modes, as noted above. 

For the Ritz approximation using 15 Ritz vectors, the junction displace- 
ments were assumed to have the form 
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r \ 

V, 


v j - |_1 cos(6 i ) Sln(3 1 ) cos(2p i ) sin(2e.) 


(82) 


and similarly for and § i . Thus, the number of junction coordinates was 
reduced from 27 to 15. In this problem the junction coordinates are far more 
important than the substructure normal mode coordinates, and it would have 
been better to retain all 27 junction coordinates and to delete all substruc- 
ture modal coordinates. However, the present example should serve to emphasize 
that reduction of junction coordinates, by any strategy, requires careful 
consideration of the particular problem at hand. 

For the uniform ring an analytical solution is available (ref. 12). 

There are three rigid-body modes. These may be given by: 


Mode 

1 : 

0 -j » o, 

w^ = cos 

V 

< 1 
II 

- sin 

Mode 

2: 

§ i = 

w. = sin 

V 

< 1 
II 

cos 

Mode 

3: 

w i = 0, 

^ = 1, 


§ i = 

- 1/R 


The first two are translation modes; the third represents rotation about the 
axis of the ring. Of course, any linear combination of the above rigid-body 
modes is also a rigid-body mode. 

The elastic in-plane bending modes can be expressed in terms of a trigo- 
nometric series in Bj. A simple solution is obtained if it is assumed that. 
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for flexural vibration, the ring is inextensional (ref. 12). In this case, 
the eigenvalues of the ring are given by 


El n 2 ( 1 - n 2 ) 
AR 4 (1 + n 2 


where n = 2 corresponds to the lowest flexural elastic mode. 

Table 6 gives the frequencies obtained for the uniform ring. The 
solutions based on Eq. (23) and on Eq. (67) are compared with the analytical 
solution given in Eq. (83). 


^^Method 

Elastic\^^ 
Mode N 

All 27 Junction 
Coordinates 
Eq. (23) 

15 Ritz Junction 
Coordinates 
E.q. (67) 

Analytical 
Eq. (83) 

1 

2.6558 x 10 3 

2.6558 x 10 3 

2.6541 x 10 3 

2 

2.6558 x 10 3 

2.6558 x 10 3 

2.6541 x 10 3 

3 

2.1249 x 10 4 

8.3493 x 10 5 

2.1233 x 10 4 

4 

2.1249 x 10 4 

8.3493 x 10 5 

2.1233 x 10 4 

5 

7.8145 x 10 4 

1.4456 x 10 6 

7.8061 x 10 4 

6 

7.8145 x 10 4 

1.4456 x 10 6 

7.8061 x 10 4 

7 

2.0450 x 10 5 

2.3357 x 10 6 

2.0416 x 10 5 


Table 6. 


Eigenvalues (u> ) 


of a Uniform Ring 
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From Table 6 it can be seen that the coupled system solution that 
includes all 27 junction coordinates produces accurate frequencies, while only 
the first two elastic mode frequencies are computed accurately by the Ritz 
representation of junction coordinates. This occurs, of course, since the 
exact mode shapes of higher-frequency modes are orthogonal to the Ritz vectors 
included in Eq. (82). As would be expected, the mode shapes of elastic modes 
1 and 2 of the Ritz solution are composed almost entirely of cos(23j) and sin 
sin(2g i ) terms. Other terms are negligible in these modes. 

Table 7 gives the frequencies of the non-symmetric ring. A solution of 
Eq. (23) containing all 27 junction coordinates is compared with results 
obtained using the Ritz approximation of Eq. (67) with 15 Ritz vectors (as in 
Eq. (82)) and 21 Ritz vectors. The latter solution included cos(3e.) and 
sin(3e i ) terms for , w i , and e. . 

Note first, from Table 7, that the coupled solution in column one 
reflects a slight raising of frequencies due to the stiffening of the ring 
and that the frequencies no longer occur in pairs. Secondly, note that neither 
of the Ritz approximations of junction coordinates is any good even though, 
in column three, the number of Ritz coordinates is almost 80% of the total 
number of physical junction coordinates. Perhaps for the ring it would be 
possible to choose more "efficient" Ritz vectors than those represented in 

Eq. (82). However, it can be argued that these are certainly logical for the 
present problem. 
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Method 

Elastic 

Mode 

All 27 Junction 
Coordinates 
Eq. (23) 

15 Ritz Vectors 
Eq. (67) 

21 Ritz Vectors 
Eq. (67) 

1 

3.3127 x 10 3 

5.6348 x 10 3 

4.1139 x 10 3 

2 

3.8962 x 10 3 

6.0030 x 10 3 

5.4910 x 10 3 

3 

2.8204 x 10 4 

1.2210 x 10 6 

3.7653 x 10 4 

4 

3.4333 x 10 4 

1.2488 x 10 6 

4.1730 x 10 4 

5 

1.0854 x 10 5 

2.2502 x 10 6 

6.1599 x 10 5 

6 

1.2269 x 10 5 

2.2889 x 10 6 

6.9169 x 10 5 

7 

2.8513 x 10 5 

3.5847 x 10 6 

1.2342 x 10 6 


Table 7. Eigenvalues (u> 2 ) of Non -axi symmetric Ring 

Equation (81) was applied to determine whether the addition of a 
cos (33^ ) contribution to the radial displacements, , would produce a signi- 
ficant change in the frequency of the third elastic mode considering the solu- 
tion with 15 Ritz vectors to be the base problem. A value of (Aw /u> 0 ) of order 
10 4 was obtained, indicating that this contribution would be very significant. 

A numerical value this high obviously violates any "small perturbation" assump- 
tion on Ao^/to 2 , but it does signify that the base problem solution for this 
mode would be too inaccurate to be useful. As a matter of fact, the solution 
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I > 


based on 21 Ritz vectors shows the cos(3g.) and si n (3e i ) contributions to 
.j , , and e.j to be dominant in the third elastic mode. 
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IMPROVEMENT OF FREE- INTERFACE METHODS 


In this part of the report attention is directed toward ways of improving 
the accuracy of free-interface methods. Example problems are solved and the 
results are compared with those of a fixed-interface method and a commonly 
used free-interface method. 


Substructure Equations 

As in Eq. (3) the substructure equation of motion may be written in 
partitioned form as 


m. . 

m. . 


* i 

X . 


k.. 

k. . 


' > 

x. 


0 

n 

1J 


1 



ij 


1 





- 


► + 



4 


► “ 4 

► 

m . . 

m . . 


X . 


k. . 

k.. 


X. 


f- 

L j 1 

JJ _ 


1 J J 


L j 1 

JJ _ 


w ^ , 


L J , 


(84) 


The free-interface modes of the substructure are obtained by setting f. = 0 
in Eq. (84) and solving the resulting equation 

mx + kx = 0 (85) 


The modes, <|> , may be normalized so that 


<f> md> = 1 

v n Y n 


(86) 


^n k< ^n = X n 


(87) 
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The modes are assembled to form the modal matrix, $, which defines substruc- 
ture generalized coordinates through the equation 

x = $ p (88) 

The modal matrix may be partitioned as follows: 


$ 


$ 


$ 


ik $ ia 
jk $ ja 


(89) 


where "k" stands for kept and "a" stands for approximated. That is, the full 
modal information of the "k" modes will be kept, but the "a" modes will be 
approximated. 

When Eq. (88) is substituted into Eq. (84), there results a set of 
uncoupled equations of the form 

P„ + P„ * 9 n (t) (90) 

where 

9„(t) * *] n fj (91) 

where <f>j n is the partition of the nth substructure mode associated with junc- 
tion coordinates (since, for free vibration the only external forcing of a 
substructure is through its interface). 

Consider now harmonic motion, as would be present in free vibration. 

Let a bar over a symbol indicate the magnitude, e.g.. 
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( 92 ) 


P n (t) = P n cos (tot) 
Equations (90) and (92) lead to 


P 


n 



(93) 


The physical displacements under harmonic motion can be obtained by combining 
Eqs. (88) and (93) to obtain 


r \ 



(94) 


where <f> n is the nth mode shape and N is the total number of modes of the 
substructure. 

Equation (94) may be written in a form which suggests the basis for the 
approximation methods of MacNeal (ref. 7) and Rubin (ref. 6), namely 


f \ 


r \ 


N 

* n + Li 


,2 2 

2 2 

x n ’ “ J 

n=S+l 

* n - a, 

V. / 


(95) 


If the second series in Eq. (95) corresponds to modes for which A n » «, 
this series may be approximated by the static effects only, and Eq. (95) 
becomes 


x 


( 


r N 


N 

^n + £ 

n n=S+l 

9n 

A 2 - \ 

X n 

n J 


l n J 


(96) 
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This corresponds to the concept of "residual flexibility" discussed by 
Klosterman (ref. 8 ). The point stressed by Rubin (ref. 6 ) is that the modes 
and frequencies indicated in the second series of Eq. (96) do not need to be 
obtained explicitly. This will be explained in the following paragraphs. 

The discussion will be restricted to harmonic motion of undamped structures 
since this is the situation encountered in the determination of modes and 
frequencies of a system. 

If the "dynamics" were neglected in Eq. (84), a "pseudo-static response, 
x^, could be obtained by solving the equation 

k x f = f ( 97 ) 

If the substructure has rigid-body freedom, k will be singular and the solu- 
tion of Eq. (97) for x f requires special consideration of this fact. Thus, 
the first and second-order approximation methods will be described first for 
constrained substructures, for which k is non-singular, and then the approxi- 
mation methods will be developed for free substructures. 

Constrained substructures - first-order method . - For harmonic forcing, 
f = f cos( w t), Eq. (97) gives 

k x ( } } = f (98) 

or 

x ( } } = G f (99) 

where G is the flexibility matrix and the superscript (1) identifies this as 
the first-order response. Since, for free vibration of a substructure. 


51 



external forces act only at junction coordinates, only the columns of G 
associated with junction coordinates are required. 

In Eq. (96) the second series represents the pseudo-static response of 
the higher-frequency modes (n >_S + 1), whereas x^ in Eq. (99) contains 
contributions from all modes. It is assumed that the modes to be kept, 
together with their natural frequencies, are available. Let Eq. (89) be 
wri tten 



$ = 

[ 

*k *a 1 

(100) 

then 

the pseudo-static response, based on Eq. (88) is 



x (1) 
x f 

= 

$ k d^P + $ d^ 

k k a a 

(101) 

or 


X 1 
— b — ' 

= 

W ♦ s*i> 

(102) 

where 


x (1) 
x k 

= 


(103) 

and 



= 

» d* 1 ) 
a a 

(104) 


Since the modes in $> are not, in general, available, it is necessary 

a 

to use Eqs. (99) and (102) to obtain x^ . To that end, let Eq. (101) be 

substituted into Eq. (98) and let the resulting equation be premultiplied 

u T . . 
by « , giving 
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[k] O k * a ] 


R’1 


r ti 

< 

* — 


d»> 


J 

i a j 


a 


{ f } 


0 05) 


Since all of the modes are mutually orthogonal with respect to k. 


*1 k «a 


« a k » k = 0 


(106) 


Also, from Eq. (87) 


*1 k 4 k 


= A,. 


(107) 


and 


K k * = A 

da a 


(108) 


where A^ and A a are diagonal matrices composed of the respective values of 
2 

. Thus, from Eq. (105) is obtained a set of uncoupled equations expressed 
collectively by 


\ ■ •! ? 


(109) 


Thus, 


.0) A-l J f 
1 k - A k *k f 


( 110 ) 


Equations (103) and (110) may be combined to give 


x (1) = G (1) f 

x k b k T 


( 111 ) 
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where 


G 0> - * k A - 1 (112) 

From Eqs. (99), (102), and (111) is obtained an expression for the contribution 
of the approximated modes to the pseudo-static response, namely 

x (1) = G ( ^ f (113) 

a a 

where 

G (1) = G - G ( J } (114) 

cl K 

Equation (113) expresses what Rubin calls the first-order contribution of the 
residual, or approximated, modes. 

The first-order method, which will also be referred to as the MacNeal 
method, may be summarized by the following equations: 

<-“ 2 1 + v p k ■ *!k (,,5) 

which is obtained from Eq. (90) for harmonic motion. Equation (96) may be 
written in matrix form as 

x (1) = \ P k + G ( a } f ( ]16 ) 

This equation will be needed when substructures are coupled to form the com- 
plete structure. In that case, the physical displacements of junction 
coordinates are required. These are given by the appropriate partition of 
Eq. (116), namely 
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(117) 


x ( P = P,, + G ( ]] f. 


jk 


JJ J 


where G^l is taken from the appropriate rows and columns of G^. 

J J cl 

Eqs. (115) and (117) are the relevant substructure equations. 


Thus , 


Constrained substructures - second-order method . - For harmonic motion, 
Eq. (84) can be written 


kx = f + w mx 


(118) 


Rubin (ref. 6) obtains the second approximation by putting the pseudo-static 
first approximation, x^^ , on the right had side of Eq. (118) and solving 

for the x on the left hand side. Thus, the pseudo-static second approximation 
is given by 


x ( ^ } = G(f + a, 2 m x ( p) 


As in Eqs. (101 ) and (102) 


j(2) 

X f 


or 


:( 2 ) 


4 >. d^ + $ d^ 
k k a a 


;< 2 k > - ;<*> 


Equations (99) and (119) can be combined to give 


3(2) = 


G (I + a) m G) f 


(119) 


( 120 ) 

( 121 ) 


( 122 ) 


By analogy with the previous derivation of Eq. (Ill) 

x (2 k } = ) 


(123) 
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= G^p (I + w 2 m G) f 


Thus, the second-order approximation of the residual (approximated) modes is 


given by 


~( 2 ) -( 2 ) 
x f x k 


G ( 2 } f 

a 


where 


• ( 2 ) = 


( G - G'^ ) ( I + a) m G ) 


This can be simplified as follows: 


( G - G ( p ) ( I + J m G ) = 


G^P [ I + J m ( G ( P + G ( P ) ] = 


Gp) + . 2 GiD.filJ) 

a a a a k 


G^P m G^P = $ Aj $T m A, ^ $>T 

a k aaakkk 


from Eq. (112). Due to orthogonality, the above vanishes. Hence, 


G<P = G<P ( I + u) 2 m G^) ) 


(127) 
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or 


where 



a 




(128) 

(129) 


The second-order approximation may be summarized by the following 
equations: 


( - O ) 2 1 + A k ) p k = *J h f- (130) 

which is the same as Eq. (115), is the primary equation for determining p k 
To effect coupling, the physical displacement equation is required, namely 

* (2) = *k ^ * G< a > f 03') 


The corresponding junction point displacements are given by 


p, + G* 2 ? f. 


x (2) ■ ■■ 

Jk k jj j 


(132) 


where G^ 2 j is the appropriate partition of G^. 


Unconstrained substructures . - If a substructure possesses one or more 
rigid-body degrees of freedom it will be referred to as an unconstrained sub- 
structure. Since k is singular, the flexibility matrix, G, used in the 
preceding analysis does not exist. On the other hand, since aerospace struc- 
tures by nature possess rigid-body freedom, it is important that methods be 
developed to analyze such structures by substructure methods similar to the 
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foregoing ones. Rubin's paper (ref. 6) was devoted to this case. 

For a substructure possessing rigid-body freedom, let the displacement 
be given by 

x = x r + x f (133) 

where is rigid-body motion and x^* is (relative) flexible motion. Also, 
let the modes of the structure be separated into rigid-body and flexible 
modes such that 

X = $ r P r + * f P f 0 34 ) 

In the notation of preceding sections all of the rigid-body modes, i.e., <t> r , 
would be kept as would a limited number of flexible-body modes. The remaining 
flexible-body modes would be approximated. Since x^ cannot be obtained by 
inverting k, as was done in going from Eq. (98) to Eq. (99), it is necessary 
to "remove" the rigid-body motion. 

Let Eq. (134) be substituted into Eq. (84) and the resulting equation 
be pre-mul ti pi i ed by to give the uncoupled equations 

p = f (135) 

^r r 

and 

P f + A f P f = f 0 36 ) 

where the orthogonality of modes is employed and it is assumed that modes are 
normalized so that 

$ r m $ r = * (137) 
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and 


<2>| m <iy 


I 


0 38) 


From Eq. (135) it is seen that, as far as rigid-body motion is concerned, 
the external forces are reacted by rigid-body inertia forces. Hence, if 
these rigid-body inertia forces are subtracted from the applied force, f, 
there will be no excitation of rigid-body modes. The physical displacement 
due to rigid-body modes is 

x r = V p r 039) 

so the (reversed) inertia force vector is given by 

f i = " m x r = " m f (140) 

Thus, the net force producing flexible-body motion is given by 

f f = f + f i = A f (141) 

where 

A = I - m $ r (142) 

To determine the flexible-body motion resulting from f^ it is necessary 
to impose temporary constraints on the substructure so that rigid-body motion 
is prevented. Let G c be the flexibility matrix of the structure which has 
been subjected to arbitrary statically-determinate constraints. G c may be 
obtained by taking k (which is singular), deleting rows and columns correspond- 
ing to the constrained degrees of freedom, inverting to get a flexibility 
matrix, and expanding the flexibility matrix by inserting zeros in rows and 
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columns corresponding to the constraints. The resulting augmented flexibility 
matrix is called G c . 

The flexible-body displacement relative to the imposed constraints is 
given by 


x<;’ - g c tf ♦ v 


(143) 


= G A f 
c 


To obtain the flexible-body displacement, x^ , rigid-body motion is removed 
from x™ by making x^ orthogonal to all rigid-body modes, i.e., 

$ T m X <1> = 0 (144) 

r f 


where 

Thus, 

or 



= x<>> + * 


C 

r r 


L<’> - * 

c r r c 


,(D 


A T x<>> 

c 


(145) 

(146) 

(147) 


From Eqs. (143) and (147) 

X W = (A T G c A) f (148) 

Equation (148) and Eq. (99) serve the same purpose for free substructures and 
constrained substructures respectively. Thus, let the symbol G be defined, 
for free substructures, by 
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G 


0 49) 


= A T G A 
c 

With this definition of G for free substructures, the remainder of the deriva- 
tion of first and second-order methods proceeds as before, beginning with 
Eq. (99). 


Coupling of Substructures 

In the preceding sections the substructure equations have been obtained 
for both constrained and free substructures. In the present section equations 
will be developed for coupling substructures. Figure 8 shows two substruc- 
tures, together with the junction coordinates at which compatibility must 
be enforced. Since there are differences in the substructure equations of 
the first-order and second-order methods, these will be treated separately. 

Coupling for the first-order approximation . - Equations (115) and (117) 
are the pertinent substructure equations for the first-order approximation. 

For substructure A they can be written 

x( jA = *jkA PkA + e( jjA fjA ( ] 50) 

and 

( - » 2 1 + A W ) p kA = f. A (151) 

For substructure B 

x( jB = $ jkB PkB + e( jjB *jB ( ] 52) 
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and 


1 + A kB ^ P kB 


J f 
jkB T jl 


(153) 


The junction coordinates of A and B must coincide. However, the number of 
"kept" modes of A and "kept" modes of B may be different. 


The coupling equations are 


(154) 


f jA = - f jB 


(155) 


where it is assumed that any required coordinate transformations have been made 
so that the junction coordinates of A coincide with the junction coordinates 
of B as given by the compatibility equation, Eq. (154). Equations (150), 

(152) and (155) may be substituted into Eq. (154) to give 


F jA ~ ^ G jjA + G ^j jB ) ( $ jkB P kB ' *jkA p kA ^ 


(156) 


Equation (156) may now be substituted into Eq. (151) to give 


( ' “ 1 + A kA ^ p kA = $ jkA ^ G ^jjA + G ^j jB ^ ^ *jkB P kB " $ jkA p kA ^ 


(157) 


Similarly, Eqs. (155) and (156) may be substituted into Eq. (153) to give 


kB ’ ^kB 


' $ jkB ^ G ^jjA + G ^jjB * ^jkB p kB “ $ jkA p kA ) 


(158) 
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Equations (157) and (158) may then be combined to form the equation 


1 

CD 

1 


r ' 

PkA 

2 

f N 

PkA 


- 


* - O) - 


_ k ba k bb 


1 P kB j 

! 

. ^ kB . 


where 


where 


K AA = A kA + $ jkA k(1) *j k A 
K AB = " $ jkA k(1) $ jkB 


K BA = " $ jkB k(1) *jkA 


K BB A kB + $ jkB k ^ °jkB 



( G 


( 1 ) 

jjA 


+ G 


( 1 ) 

jjB 



( 159 ) 


( 160 ) 


(161) 


Equation (159) thus forms the primary eigenvalue problem for determining the 
modes and frequencies of the coupled structure. It may be noted that its 
size is equal to the sum of the kept substructure modes in substructures 
A and B. 


Coupling for the second-order approximation . - The coupling equations 
for the second-order approximation can be obtained from Eqs. (157) and (158) 
by replacing G^jj by G^j for the substructure. The complication arises due 
to the fact that G^l is a submatrix of G^, which is given by Eq. (128) as 
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(128) 


Thus, 


G< 2 > = + „2 H 0) 

a a a 


G^l. 

JjA 


= G 


( !1. +U, 2 h ( 1! 


jjA 


jj A 


(162) 


and 


r (2) 

G jjB 


g<!> + » 2 h<!> 

JjB JjB 


(163) 


Let 


r (l) + r (D 

G jjA + G JjB 


= G 


( 1 ) 

jj 


(164) 


and 


u(l) + u(l) 

H jjA H jjB 


H ( l? 

JJ 


(165) 


Then, 


g ( 2 > + g< 2 > = g<!! +u 2 h<!) 

JjA J JB jj jj 


(166) 


From Eqs. (157). and (158) the system equations for the second- 
order method may be deduced to be 


A + k ( ^ ) a 

\A jkA K jkA 

_ $ |<(2) $ 
jkA K ®jkB 

1 

r n 

PkA 

2 

’ = to 4 

r \ 

PkA 

' ‘JkB k(2) ‘jkA 

_ 

(21 T 

A kB + *jkB k $ jkB 


Kb 

\ > 


PkB 

l J 


(167) 
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where 


.( 2 ) 


( G 


( 2 ) 

jjA 


+ G 


( 2 ) 

jjB 


-1 


( G 


(1) 

jj 




(168) 


2 

Since, from Eqs. (166), oo appears on both the left and right-hand sides 
of Eq. (167) this system equation cannot be solved conveniently in its present 
form. Two methods may be used to circumvent this difficulty: (1) an iteration 

method may be used, with the u on the left-hand side of Eq. (167) being taken, 
in the first iteration, as the u> obtained by using the first-order method, 
or (2) a Maclaurin series expansion may be employed to obtain an approximation 
to the inverse matrix required on the left-hand side of the system equation. 

The Maclaurin series method will be described more fully now. Let 


' 2 > - ( g<!> ♦ „, 2 «<!> r 


JJ 


« H' ) 


JJ 


( G<]> [ I ♦ » 2 (g'J])- 1 H<]> ] )■’ 


[ i + w 2 icd) h ( ! 1 ]'' kd) 

J J 


(169) 


From Eq . (169) it can be seen that the second term in the square bracket dis- 
tinguishes the second-order method from the first-order method. If the correc- 
tions are "small" k^ can be approximated by 


.( 2 ) 


( I 


kd)„d). 

J J 


) k 


( 1 ) 


(170) 
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From Eqs. (167) and (170) the following second-order system equation is 



( PkB 


where 


M M ' 1 + 4 jkA I "* 1) ^jkA 


"JkA mll) { jkB 


- *JkB ra<1) *jkA 


H BB ' 1 + *jkB *jkB 


where 


,0) . k (l) „0) k (D 
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Examples 


Several example problems were solved by the free-interface methods 
described in the preceding sections. The purpose of these examples is to 
illustrate the results obtainable by the first-order and second-order methods 
and to compare these results with those obtained by using a fixed-interface 
method (method of Craig and Bampton) and another free-interface method 
(Hou's method) . 

Beam examples - first and second-order methods versus Hou method . - 
Frequencies of a two-substructure uniform beam were obtained for three sepa- 
rate sets of boundary conditions: clamped-clamped , clamped-free, and free-free. 
Thus, the theory for both constrained and free substructures was employed in 
these examples. Figure 9 shows the configuration of the two substructures for 
each of the three support conditions. Substructure A consists of four elements, 
while substructure B consists of three elements. Consistent mass matrices 
were employed for element mass matrices. 

In Table 8 are shown the system frequencies obtained for the clamped- 
clamped beam by Hou's method and by the first and second-order methods. To 
obtain each frequency by the second-order method, the iteration procedure 
was employed, wherein the corresponding frequency obtained by the first-order 
method was used in determining the of Eq. (168). 

Note that the first-order (MacNeal) method is much more accurate than 
Hou's method, and that further improvement is obtained by use of the second- 
order (Rubin) method. On the other hand, increasing the number of substruc- 
ture modes while using the first-order method also produces significant 
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improvement. 

The "exact" values in Table 8 were obtained by retaining all original 
physical coordinates. 

Table 9 shows a more limited set of results for the clamped-free beam. 
The accuracy of the first-order method is not as good as it was for the 
clamped-clamped beam; but the first-order method is still far superior to 
Hou's method. 

Table 10 gives results obtained for the free-free beam. Note that the 
results obtained by the first-order method are virtually identical to those 
•obtained for the clamped-clamped beam. Of course, the free-free beam has two 
zero-f reguency rigid-body modes which are not listed in the table. Again 
Hou's method does not produce acceptable accuracy. 




— J 




Elastic 

Mode 

Number 

II 

II 

CO 

z: 

N a = 4 

» N B " 4 

Exact 

1st order 
(MacNeal ) 

% 

error 

Hou 

% 

error 

1 

.0000052 

0.00 

.0000053 

1.92 

.0000052 

2 

.0002069 

2.27 

.0002624 

29.71 

.0002023 

3 

.0016227 

2.15 

.0017612 

10.86 

.0015886 

4 

.0061338 

0.01 

.0072102 

17.56 

.0061330 

5 

.0177475 

4.73 

.0228454 

34.81 

.0169463 

6 

.0385308 

0.21 

.0418527 

8.85 

.0384504 

7 

.0828725 

10.58 

— 

-- 

.0749409 

8 

.5902060 

X 

— 

-- 

.1602931 


Table 9. Frequencies of a Two-Substructure Clamped-Free Beam 


Elastic 

Mode 

Number 

N A = 5 ’ 

n b m 

N a =5 

N b M 

Exact 

1st order 
(MacNeal) 

% 

error 

Hou 

0/ 

/o 

error 

1 

.0002086 

0.05 

.0002807 

34.63 

.0002085 

2 

.0015875 

0.02 

.0017529 

10.44 

.0015872 

3 

.0061381 

0.14 

.0072018 

17.49 

.0061297 

4 

.0170057 

0.55 

.0228352 

35.02 

.0169130 

5 

.0385203 

0.69 

.0415545 

8.62 

.0382562 

6 

.0783702 

5.52 

— 

— 

.0742702 

7 

.3928161 

X 

— 

— 

.1573651 


Table 10. Frequencies of a Two-Substructure Free-Free Beam 
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Truss examples - MacNeal method versus Craig-Bampton method . - The 
Benfield truss of Fig. 3 was used in a study comparing the results obtainable 
by the first-order method with results obtained by use of the Craig-Bampton 
method. The results are shown in Table 11. The exact frequencies are the 
same as those given in Table 1, and were obtained using all of the original 
physical coordinates of the coupled system. 

Since the Craig-Bampton method employs junction coordinates as well as 
substructure normal mode coordinates, the first-order method and Craig- 
Bampton method are compared on the basis of total number of system modes 
(N tot = 16). Since no convergence studies were carried out, no strong con- 
clusions can be reached about the comparison of MacNeal 's method with the 
Craig-Bampton method, but Table 11 suggests that the MacNeal method produces 
results of accuracy comparable to those produced by the Craig-Bampton method. 
The Craig-Bampton method produces the best results for the low-frequency 
modes. The MacNeal method is not quite as accurate for these modes, but 
does better for the modes in the middle range of frequencies. 



Elastic 

Mode 

Number 


n a = 

n b = s 

1st order 
(MacNeal) 

% 

error 


0.25 

0.26 

1.60 

0.36 

0.64 

0.97 

2.01 

2.76 

0.28 

3.54 

33.04 

75.15 


n a s n b ■ 

5, Nj = 6 

Exact 

HI 

% 

error 


0.004391 

0.00 

0.004391 

0.018331 

0.01 

0.018329 

0.030527 

0.06 

0.030508 

0.041664 

0.18 

0.041589 

0.068734 

0.22 

0.068585 

0.098938 

0.17 

0.098770 

0.126051 

10.61 

0.113959 

0.144412 

18.02 

0.122362 

0.161327 

7.05 

0.150701 

0.247411 

58.82 

0.155147 

0.361008 

X 

0.174052 

0.506162 

X 

0.200979 

1.030503 

X 

0.215105 


0.004402 

0.018377 

0.030997 

0.041740 

0.069022 

0.099727 

0.116246 

0.125738 

0.151118 

0.160644 

0.231554 

0.352010 

1.524661 


Table 11. Frequencies of Two-Substructure Benfield Truss by 
MacNeal Method and by Craig-Bampton Method 
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Grid examples - MacNeal method versus Craiq-Bampton method for 
cantilever and free-free grid structures . - The grid structure shown in 
Fig. 5 was used for further comparison of the MacNeal method and the Craig- 
Bampton method. A similar free-free grid was also analyzed. In these studies 
a consistent mass matrix was employed. 

Table 12 shows the results obtained for the first twenty modes of the 
cantilever grid. The column labelled "exact" is for the 72 d.o.f. original 
structure. For smaller numbers of total modes, the MacNeal method produces 
considerably better results than the Craig-Bampton method. 

Figure 10 shows the number of modes accurate to within 0.5% as computed 
by each of the two methods. (The results shown in Fig. 10 include all the 
"0.5% modes" among the total number of modes calculated, i.e., 20, 24, 28, 

32, even though only twenty modes are tabulated in Table 12.) 

Table 13 shows the results obtained for the first twenty frequencies 
of a free-free grid structure as obtained by the MacNeal method and by the 
Craig-Bampton method. No exact values are available because of limitations 
on the size problem which could be handled by the eigensolver being used. 

It is noted that the frequencies are quite close, for the most part, with 
the maximum difference being only about 2.5%. 
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EXACT 

N tot ■ 20 
V 11 , Ng = 9 

N tot ! 

V 5 . Ng = 3 

20 

V 12 

MacNeal 

% 

error 

Craig-Bampton 

% 

error 

.00006210 

.00006210 

0.00 

.00006211 

0.01 

.00025109 

.00025110 

0.00 

.00025115 

0.02 

.00163979 

.00164073 

0.05 

.00164660 

0.41 

.00238797 

.00238799 

0.00 

.00239788 

0.41 

.00394372 

.00394463 

0.02 

.00395306 

0.23 

.00819127 

.00820367 

0.15 

.00871273 

6.36 

.00862821 

.00865644 

0.32 

.00881859 

2.20 

.01806006 

.01808477 

0.13 

.01857390 

2.84 

.01991494 

.01998109 

0.33 

.02188961 

9.91 

.02254916 

.02260244 

0.23 

.02271096 

0.71 

.02744331 

.02752603 

0.30 

.03864739 

40.82 

.03262871 

.03298812 

1.10 

.06218932 

90.59 

.04493968 

.05362282 

19.32 

.08120924 

X 

.05150864 

.05581855 

8.36 

.12613476 

X 

.06175304 

.06373244 

3.20 

.20728453 

X 

.06664190 

.06884813 

3.31 

.36264343 

X 

.06963443 

.17950904 

X 

.53033012 

X 

.08534049 

.25447977 

X 

1.05181024 

X 

.08680764 

.69330843 

X 

2.45246230 

X 

.08914933 

1.32773903 

X 

5.23913299 

X 


Table 12. Frequencies of Cantilever Grid Structure 
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Elastic 

Mode 

N tot' 
N a =14 , 

24 

n b = io 



28 

N b = 12 

Number 

MacNeal 

% 

error 

Craig-Bampton 

% 

error 

MacNeal 

l 

error 

1 

.00006210 

0.00 

.00006210 

0.00 

.00006210 

0.00 

2 

.00025110 

0.00 

.00025109 

0.00 

.00025109 

0.00 

3 

.00164035 

0.03 

.00164019 

0.02 

.00164025 

0.02 

4 

.00238799 

0.00 

.00239012 

0.09 

.00238797 

0.00 

5 

.00394406 

0.00 

.00394433 

0.01 

.00394400 

0.00 

6 

.00820253 

0.13 

.00819496 

0.04 

.00819978 

0.10 

7 

.00863702 

0.10 

.00867873 

0.58 

.00863315 

0.05 

8 

.01807972 

0.10 

.01825831 

1.09 

.01806091 

0.00 

9 

.01996598 

0.25 

.01997539 

0.30 

.01993582 

0.10 

10 

.02257782 

0.12 

.02267290 

0.54 

.02255448 

0.02 

11 

.02750739 

0.23 

.02750874 

0.23 

.02750444 

0.22 

12 

.03289693 

0.82 

.03320131 

1.75 

.03266930 

0.12 

13 

.04507813 

0.30 

.05276010 

17.40 

.04505674 

0.26 

14 

.05574645 

8.22 

.05686502 

10.39 

.05158030 

0.13 

15 

.06250395 

1 .21 

.06696988 

8.44 

.06177593 

0.03 

16 

.06861338 

2.95 

.08138276 

22.11 

.06858598 

2.91 

17 

.07109523 

2.09 

.27636744 

X 

.06969506 

0.08 

18 

.08733079 

2.33 

.35125316 

X 

.08557656 

0.27 

19 

.09138064 

5.26 

.39979778 

X 

.08701777 

0.24 

20 

.10946900 

22.79 

.53601805 

X 

.09134459 

2.46 


Table 12 (cont.). Frequencies of Cantilever Grid Structure 
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Elastic 

Mode 

NS ' ' *v'- 1< J ; " 

N tot ■ 3 
N A =18 , N 

ET 14 

N tot * 32 
N A = 12 ,N B = 8 , 

N. = 12 

J 

Number 

Craig-Bampton 

% 

error 

MacNeal 

% 

error 

Craig-Bampton 

% 

error 

1 

.00006210 

0.00 

.00006210 

0.00 

.00006210 

0.00 

2 

.00025109 

0.00 

.00025109 

0.00 

.00025109 

0.00 

3 

.00163996 

0.01 

.00163998 

0.01 

.00163989 

0.00 

4 

.00238879 

0.03 

.00238797 

0.00 

.002.38814 

0.00 

5 

.00394429 

0.01 

.00394398 

0.00 

.00394381 

0.00 

6 

.00819481 

0.04 

.00819423 

0.03 

.00819336 

0.02 

7 

.00863617 

0.09 

.00863294 

0.05 

.00863058 

0.02 

8 

.01822469 

0.91 

.01806079 

0.00 

.01806936 

0.05 

9 

.01995132 

0.18 

.01993480 

0.09 

.01992609 

0.05 

10 

.02267161 

0.54 

.02255635 

0.02 

.02255262 

0.01 

11 

.02748344 

0.14 

.2746659 

0.08 

.02746586 

0.08 

12 

.03313560 

1.55 

.03266878 

0.12 

.03265385 

0.07 

13 

.04505185 

0.24 

.06505075 

0.24 

.04498599 

0.10 

14 

.05672962 

10.13 

.05157360 

0.12 

.05157714 

0.13 

15 

.06272967 

1.58 

.06177415 

0.03 

.06186679 

0.18 

16 

.06695281 

0.46 

.06735587 

1.07 

.06694378 

0.45 

17 

.07023630 

0.86 

.06967776 

0.06 

.06967930 

0.06 

18 

.08792538 

3.02 

.08553947 

0.23 

.08547908 

0.16 

19 

.08987727 

3.53 

.08699608 

0.21 

.08694733 

0.16 

20 

.12568251 

40.97 

.09114901 

2.24 

.08962024 

0.52 


Table 12 (cont.). Frequencies of Cantilever Grid Structure 
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Elastic 

Mode 

Number 

MacNeal 

N a = 20 , N b = 15 
N tot ■ 35 

Craig-Bampton 
N A = 12 , n b = 8 , n. = 15 

N tot - 35 

1 

.00045914 

.00045913 

2 

.00091210 

.00091192 

3 

.00197905 

.00197929 

4 

.00286973 

.00286953 

5 

.00399206 

.00399188 

6 

.00680185 

.00679642 

7 

.01032129 

.01032054 

8 

.01156372 

.01156373 

9 

.01462129 

.01472953 

10 

.01979609 

.01986503 

11 

.02500327 

.02500680 

12 

.02542485 

.02538976 

13 

.03188867 

.03237293 

14 

.03479776 

.03476012 

15 

.05306692 

.05435942 

16 

.05514919 

.05609568 

17 

.05588742 

.06071017 

18 

.06206941 

.06240902 

19 

.06640768 

.06484146 

20 

.08028079 

.07861373 


Table 13. Frequencies for Two-Substructure Free-Free Grid 




Testing Requirements 


The free-interface methods described in this report, in particular the 
first-order method, are compatible with dynamic test procedures, as was 
pointed out by Rubin (ref. 6). Although considerable work needs to be done 
in exploring the relationship of these methods to dynamic test procedures, 
several comments may be made on the basis of the equations presented in this 
report. 

The first-order method requires the substructure information indicated 
in Eqs. (115) and (117): 

(D A k - square of natural frequencies of substructure free- 

interface modes to be kept 

(2) <5^ - normalized (with respect to the mass matrix) free-inter- 

face mode shapes of modes to be kept 

(3) gH) - residual flexibility at junction coordinates 

0 J 

Items (1) and (2) are standard items which are determined in mode 

survey tests. The only complicating factor is that the modes are to be 

normalized with respect to the mass matrix. 

The determination of the residual flexibility at junction coordinates 

appears to be relatively straightforward for constrained substructures. 

However, since 1 is a submatrix of GpP which is determined as the differ- 
j J a 

ence of two flexibility matrices, Eq. (114), the sensitivity of this matrix 
to measurement errors should be thoroughly explored. 
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For unconstrained substructures the determination of the residual 
flexibility matrix requires a constrained flexibility matrix, G c , and a trans- 
formation matrix A, which, by Eq. (142), involves the mass matrix and rigid- 
body modes. Equations (149) defines the flexibility matrix modified to 
account for rigid-body modes. Considerable study is needed to determine: 

(a) whether Eq. (149) provides the best way to treat rigid-body modes, and 

(b) whether Eq. (149) is sensitive to measurement errors. 
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CONCLUSIONS AND RECOMMENDATIONS 


The studies of fixed-interface methods indicate that: 

(1) Reduction of the number of junction coordinates is feasible since 
this produces results of accuracy comparable to that obtained when 
substructure normal modes are reduced. 

(2) Reduction of junction coordinates may be based on: 

(a) Ritz representation of junction coordinates using Ritz 
vectors provided by the analyst. 

(b) Guyan reduction of junction coordinates using a static 
reduction of coordinates, with the retained coordinates 
selected by the analyst. 

(c) Modal reduction of junction coordinates, using "junction 
modes" as Ritz vectors and using frequency ratio or modal 
participation factor to determine the coordinates to be 
retained. 

(3) Although the choice of method for reducing the number of junction 
coordinates is dependent, to some extent , on the problem to be 
solved, modal reduction has two advantages not offered by the 
other two methods : 

(a) The calculation of modal participation factors is very simple. 

(b) The final equations are of a form that could lead to a very 
efficient eigensolution algorithm if advantage were taken of 
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the fact that the generalized stiffness matrix is diagonal. 

A disadvantage is that a preliminary eigensolution for junction 
modes is required. 

(4) Error analysis, i.e., determination of the frequency "error" intro- 
duced by omission of a particular Ritz vector (mode), may be con- 
ducted for any of the three proposed methods. The calculations are 
somewhat tedious except for the modal reduction method. 

(5) Before investing further effort in implementing any of the above- 
named methods for reduction of junction coordinates, an attempt 
should be made to create an efficient computer program for imple- 
menting the original Craig-Bampton equation, Eq. (23). It should 
be possible to create an efficient power iteration program, since 
K^. is diagonal and Kjj will be narrowly banded or, at worst, banded 
with only a few submatrices outside of the band. If such a program 
were available, it would be a straightforward matter to provide an 
option for modal reduction of junction coordinates. 

Free-interface substructure models suggested by MacNeal and Rubin have 
been studied. Equations for coupling substructures to form systems have been 
presented and a number of example problems involving two substructures have 
been solved. On the basis of these studies, it is concluded that: 

(1) Only substructure normal mode coordinates are retained in the final 
system equations of motion, yet, both first-order (MacNeal) and 
second-order (Rubin) free-interface methods produce accuracies 
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comparable to those produced by Hurty-type fixed-interface methods. 


The second-order (Rubin) method requires additional calculations, 
and further study is required to establish whether the improvement 
in results is sufficient to justify the added computational effort. 

The first-order method appears to be compatible with dynamic test 
procedures, but further study is needed, particularly with respect 
to unconstrained substructures. 

"Rubin-type" methods appear to give upper bounds to the system 
frequencies and to be rapidly convergent. Further study of error 
bounds is required, however, to confirm these points. 

Mode shape and stress errors should be investigated for "Rubin- 
type" methods and should be compared with those produced by Hurty- 
type methods. 
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FIG. I . COMPARISON OF METHODS WITH FREQUENCY ERROR OF <0.5% 



• interior coordinates (x, ) 
o junction coordinates (xj ) 


FIG. 2. SUBSTRUCTURE COORDINATES 



FIG. 3. TWO SUBSTRUCTURE BENFIELD TRUSS 
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NO. MODES ACCURATE WITHIN 0.5% 



FIG. 4. MODAL INTERFACE REDUCTION vs GUY AN 

INTERFACE REDUCTION FOR BEN FIELD TRUSS 




FIG. 5 . CANTILEVER GRID STRUCTURE 
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Case I fj= f ^ (uniform) 

Case2 f 2f / 

FIG. 7 . RING STRUCTURES 
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NO. MODES ACCURATE WITHIN 0.5% 



FIG. 9. BEAM CONFIGURATIONS 



FIG. 10. COMPARISON OF MACNEAL AND CRAIG- BAM PTON 
METHODS FOR A CANTILEVER GRID STRUCTURE 


NASA-Langley, 1977 CR-2781 
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